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© The frequency-versus-time characteristic of the 
output signal of a swept-frequency radar is linearized 
by determining departures from linearity of the 
phase-versus-time characteristic of the beat signal 
obtained between an actual return from a selected 
radar target and a sample of the current radar output 
signal, and applying oppositely-directed corrections 
to the frequency-versus-time characteristic. The rel- 
evant beat signal is selected from the output signal 
of a mixer (8) by performing spectral analysis in a 
Fast Fourier Transform (FFT) calculating circuit (15) 
and then performing an inverse FFT on a windowed 
portion of the FFT obtained. A second or higher- 
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order polynomial is then fitted to the phase-versus- 
time characteristic of the selected beat signal after it 
has been expressed as time versus phase, and 
correction terms are derived from the coefficients of 
this polynomial. The radar output signal is generated 
by a voltage-controlled oscillator (1) controlled by a 
control voltage generator (7). The generator com- 
prises a voltage source (27) and a cascade combina- 
tion of integrators (18,19,20) outputs of which are 
connected to a voltage summing circuit (13) via 
respective multipliers (32,33,34,35). The correction 
terms are applied as multiplication factors to the 
higher-order multipliers (33,34,35). 
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The present invention relates to a method of, 
and apparatus for, reducing non-linearity (if any) 
present in the frequency-versus-time characteristic 
of the output signal of a swept-frequency radar 
transmitter. 

Swept-frequency radar apparatuses which 
transmit an output signal which has a substantially 
linear frequency-versus-time characteristic during 
each frequency sweep are well-known. If the cur- 
rent output signal is mixed with a version of said 
output signal which has been transmitted and re- 
flected back by a target a beat signal results, 
because the frequency of the output signal will 
have changed during the time taken by the trans- 
mitted output signal to travel to the target and back 
again. The frequency of the beat signal arising 
from reflection by a particular target will therefore 
be representative of the range of that target if the 
velocity of the target relative to the radar apparatus 
is zero, or will be representative of a combination 
of the range and the relative velocity of the target 
otherwise, because any non-zero relative velocity 
will give rise to a Doppler frequency shift in the 
reflected signal. Thus analysis of the beat signal 
into its respective frequency components will yield 
range/relative velocity information about any target 
which has given rise to a reflection. 

It will be appreciated that a given target at 
constant range will only give rise to a constant beat 
frequency if the frequency-versus-time characteris- 
tic of the radar transmitter output signal is ac- 
curately linear. Any increase or decrease of the 
(positive or negative) rate of change of frequency 
during a given frequency sweep will give rise to an 
increase or decrease respectively in the frequency 
of each frequency component of the beat signal 
during that sweep, effectively resulting in the 
spreading of each frequency component and hence 
in a decrease in the resolution of the range/velocity 
information contained therein. 

In order to optimise the frequency sweep lin- 
earity it is known from, for example, a paper "On 
performance of a linear FM radar transmitter at 
35GHz" by P.Z. Peebles and A.H. Green at pages 
4-10 of Proc. IEEE Southeast Con. 1982 to employ 
a closed-loop feedback system to provide dynamic 
correction. To this end the output signal of the 
radar transmitter, is applied 'directly to one input of 
a mixer and also to the other input of the mixer via 
a delay line having a small delay relative to the 
repetition period of the (periodic) frequency sweeps 
of the output signal. The resulting output of the 
mixer is a beat signal which ideally (if the 
frequency-versus-time characteristic of the output 
signal is perfectly linear) has a constant frequency 
f 0 equal to the product of the delay produced by 
the delay device and the rate of change of fre- 
quency of the output signal. This beat signal is 



applied to what is in effect a frequency discrimina- 
tor tuned to f 0 and the discriminator output signal is 
combined with the (sawtooth) modulation waveform 
which produces the frequency sweeps in the output 
5 signal in such a sense as to reduce any deviation 
of the frequency of the beat signal from f 0 . Thus 
any non-linearity in the frequency-versus-time char- 
acteristic of the output signal, which non-linearity 
gives rise to a deviation in the frequency of the 
io beat signal from f c , results in a modification to the 
modulation waveform such as to reduce the non- 
linearity. The control loop in effect continuously 
maintains the rate of change of frequency of the 
transmitter output signal at a value such that the 
is frequency of the resulting beat signal is f 0 . 

A disadvantage of the known feedback system, 
at least in some potential applications, such as 
automotive radar, is the cost and/or bulk of the 
additional microwave or high-frequency compo- 
20 nents required thereby, e.g. the delay line and the 
means by which the delay line is coupled to the 
remainder of the apparatus. 

It is an object of the present invention to miti- 
gate this disadvantage. 
25 According to one aspect of the present inven- 
tion there is provided a method of reducing non- 
linearity (if any) present in the frequency-versus- 
time characteristic of the output signal of a swept- 
frequency radar transmitter, which output signal is 
30 derived from a radio-frequency source the output 
frequency of which is dependent on the value of a 
control signal applied in operation to a control 
signal input of said source, the output of a control 
signal generator being coupled to said control sig- 
35 nal input for applying said control signal thereto to 
repeatedly sweep the frequency of said output 
signal, in which method said output signal is mixed 
with a delayed version thereof to yield a beat signal 
and said characteristic is adjusted in accordance 
40 with any variation of a frequency of this beat signal 
occurring during individual frequency sweeps of 
said output signal, characterized in that said control 
signal is formed by combining a time-independent 
component and n individually weighted time-de- 
45 pendent components, where n is greater than one, 
where the first time-dependent component is de- 
rived from the output of a single resettable ana- 
logue integrator circuit which is fed from the output 
of a d.c. reference source and where the mth time- 
so dependent component (m = 2 n) is derived from 

the output of a cascade combination of m reset- 
table analogue integrator circuits which is fed from 
the output of a d.c. reference source, in that the 
delayed version of the transmitter output signal is 
55 constituted by a version of this output signal which 
has been transmitted and reflected back by a tar- 
get or targets, in that the cumulative-and- 
normalized-phase versus time characteristic of a 
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specific spectral component of said beat signal is 
determined for an individual frequency sweep of 
said transmitter output signal, in that the coeffi- 
cients of the first nth order terms of an nth order 

polynomial which fits the dependence of said con- 
trol signal on said cumulative and normalized 
phase are determined from said cumulative-and- 
normalized-phase versus time characteristic and 
the concomitant variation of said control signal with 
time, and in that the coefficients so determined are 
used to update the weighting coefficients of the 
f' r st nth said time-dependent components re- 
spectively to reduce any said non-linearity. 

It has now been recognized that it is not essen- 
tial to provide a dedicated delay line to create a 
reference signal for mixing with the current output 
signal of the radar transmitter; the transmitter out- 
put signal may itself be used for this purpose after 
transmitting it to and receiving it back from a 
suitable reflecting target. The resulting beat signal 
will, of course, no longer have a predetermined 
average frequency, (unless the same stationary tar- 
get is always employed) but changes occurring in 
the actual frequency thereof during individual fre- 
quency sweeps of the transmitter output signal will 
still be indicative of sweep non-linearity and hence 
of the adjustments required to the frequency- 
versus-time characteristic of the transmitter output 
signal. These beat signal frequency changes mani- 
fest themselves as non-linearities in the cumulative 
phase versus elapsed time characteristic of the 
beat signal during the individual frequency sweeps, 
and it has now been recognized that, if this char- 
acteristic is measured, thereby providing informa- 
tion on the basis of which the requisite adjustment 
can be carried out, the various values of cumulative 
phase in the measured characteristic can further- 
more be normalized, i.e. can each be divided by 
the value which the cumulative phase has a spe- 
cific time after the beginning of the relevant fre- 
quency sweep, for example at the end thereof, 
without loss of the required information. After nor- 
malization the same characteristic will result what- 
ever the basic beat frequency is, i.e. whatever the 
range of the corresponding target is. Once the 
cumulative and normalized phase versus elapsed 
time characteristic of the beat signal has been 
determined for an individual frequency sweep of 
the transmitter output signal, the dependence of 
the cumulative and normalized phase on the value 
of the frequency control signal for the transmitter 
radio-frequency source, and hence the dependence 
of said control signal on said cumulative and nor- 
malized phase, can be worked out from a knowl- 
edge of how this control signal varied with elapsed 
time during the same sweep. It will be evident that, 
in order to produce a beat signal of constant fre- 
quency during a subsequent frequency sweep, i.e. 



a beat signal which has a linear cumulative and 
normalized phase characteristic, the frequency 
control signal for the radio-frequency source should 
vary during that sweep in such a manner that it is 
5 the same function of normalized elapsed time as it 
has been found to be of normalized beat signal 
cumulative phase. Thus a knowledge of how the 
control signal depends on beat signal normalized 
cumulative phase during a given frequency sweep 
to provides directly a knowledge of how the control 
signal should be made to depend on normalized 
elapsed time during subsequent frequency sweeps. 

According to another aspect of the present 
invention, there is provided a radar apparatus com- 
75 prising a radio-frequency source having an output 
coupled to signal operating means for radiating a 
transitter signal and for receiving at least a return 
signal, a control signal generator having an output 
coupled to a frequency control signal input of said 
20 radio-frequency source for applying a frequency 
control signal to said input to repeatedly sweep the 
frequency of the output signal of said source, and a 
mixer having a first input coupled to the signal 
operating means for receiving a version of the 
25 output signal of said source which has been trans- 
mitted and reflected back by a target or targets, a 
second input coupled to the output of said radio- 
frequency source for receiving a sample of the 
output signal of said source as a reference signal, 
30 and an output for a beat signal between said ver- 
sion of the output signal of said source and said 
sample, said control signal generator comprising 
component combining means for combining a time- 
independent component and n individually weight- 
35 ed time-dependent components to form said fre- 
quency control signal, where n is greater than ons, 
a resettable analogue integrator circuit arrangement 
for generating said time-dependent components, 
and means for individually weighting each said 
40 time-dependent component, said component com- 
bining means having a respective input corre- 
sponding to each said time-dependent component 
and said resettable analogue integrator circuit ar- 
rangement comprising a single resettable analogue 
45 integrator circuit included between the output of a 
d.c. reference source and the combining means 
input for the first time-dependent component and a 
cascade combination of m resettable analogue in- 
tegrator circuits included~between the output of a 
50 d.c. reference source and the combining means 
input for the mth time-dependent component 

(™ - 2 n). the apparatus further including means 

for determining the cumulative-and-normalized- 
phase versus time characteristic of a specific spec- 
55 tral component of the mixer output signal for an 
individual frequency sweep of the radio-frequency, 
source, means for determining the coefficients of 
the first nth order terms of an nth order polyno- 
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mial which fits the dependence of said frequency 
control signal on said cumulative-and-normalized- 
phase from said cumulative-and-normalized-pha- 
seversus time characteristic and the concomitant 
variation of said control signal with time, and 
means for updating the weighting coefficients of 

the first nth said time-dependent components in 

accordance" with the coefficients of the first nth 

order terms respectively so determined: to reduce 
non-linearity (if any) present in the frequency 
sweeps of the output signal of said radio-frequency 
source. 

Embodiments of the invention will now be de- 
scribed, by way of example, with reference to the 
accompanying diagrammatic drawings, in which: 

Figure 1 is a block diagram of radar apparatus 

for implementing a method in accordance with 

the invention, and 

Figure 2 is a flow diagram illustrating operations 
performed by a programmed microcomputer in- 
cluded in the apparatus of Figure 1. 
In Figure 1 radar apparatus comprises a swept- 
frequency radar transmitter which includes a 
voltage-controlled r.f. oscillator 1 the output 2 of 
which is coupled to a transmit/receive aerial 3 via a 
circulator 4. The frequency control input 5 of os- 
cillator 1 is fed from the output 6 of a control 
voltage generator 7. 

A receiver portion of the radar apparatus in- 
cludes a mixer 8 the output 9 of which is coupled 
to the input 14 of a Fast Fourier Transform (FFT) 
calculating circuit 15 via an amplifier 12, a low-pass 
filter 13 and an analogue-to-digital converter 11. A 
first input 10 of mixer 8 is fed from the aerial 3 via 
circulator 4, and a second input 16 of mixer 8 is 
fed with a sample of the output signal of oscillator 
1 via a coupler 17. 

The output waveform of control voltage gener- 
ator 7 is basically a linear sawtooth. On the as- 
sumption that the output frequency of oscillator 1 is 
linearly dependent upon the control voltage fed to 
its input 5 this output frequency in consequence 
basically repeatedly sweeps from a first frequency 
F 0 to a second frequency Fi at a constant rate A. If 
the resulting signal transmitted from aerial 3 is 
reflected from an object at a range R and received 
back at the aerial 3 then, as is known, the dif- 
ference between the instantaneous, frequencies ap- 
plied to the inputs 10 and 16 of mixer 8 (the 
instantaneous frequencies of the signal received by 
aerial 3 and of the output signal of oscillator 1 
respectively) is equal to 2AR/C, where c is the 
velocity of light, possibly Doppler-shifted due to 
relative motion between the radar apparatus and 
the reflecting object. The output signal of mixer 8 
therefore has a component having this frequency 
f b . After amplification by amplifier 12, low-pass 
filtering by filter 13 and conversion to digital form 



by converter 11, this output signal is applied to 
FFT calculating circuit 15. (2N-point) FFT calculat- 
ing circuit 15, as is known, acts as a bank of N 
filters and in consequence performs a spectral ana- 
5 lysis, producing a signal on that one of its N 
outputs which corresponds to a frequency range 
which includes the frequency f b . The presence of a 
signal on this output therefore indicates the pres- 
ence of a stationary reflecting object at the range Ft 
w (or of a moving object at a range somewhat dis- 
placed from this range R). The simultaneous pres- 
ence of a plurality of reflecting objects at different 
ranges will, of course, result in signals at a cor- 
responding plurality of outputs of the circuit 15. 
75 In practice the output frequency F of oscillator 

1 will not be exactly linearly dependent upon the 
value of the control voltage applied to its input 5, 
with the result that, if the output signal of control 
signal generator 7 is a linear sawtooth, the rate of 
20 change A of this frequency will vary during the 
course of each frequency sweep. This means that 
the frequency of the component of the output sig- 
nal of mixer 8 corresponding to each reflecting 
object will itself vary during the course of each 
25 frequency sweep, resulting in signals at adjacent 
outputs of FFT-calculating circuit 15; the range 
resolution of the radar apparatus will be impaired. 

In order to improve matters in this respect the 
output waveform of control voltage generator 7 is 
30 automatically modified in such a manner as to 
compensate, at least partly, for the non-linear de- 
pendence of the output frequency of oscillator 1 on 
its control voltage, i.e. in such a manner as to 
improve the constancy of the rate of change A of 
35 the output frequency of oscillator 1 during the 
course of each frequency sweep. Each nominally 
linear ramp of the output voltage of generator 7 is 
distorted into a curve whose curvature takes the 
opposite form to the curve of output frequency of 
40 oscillator 1 against control voltage applied to input 
5. 

To enable this to be done generator 7 includes 
a cascade combination of, in the present case 
three, integrator circuits 18,19 and 20 having in- 

45 tegration signal inputs 21,22 and 23 respectively, 
and reset signal inputs 24,25 and 26 respectively. 
Input 21 of circuit 18 is fed from a d.c. voltage 
source 27, and inputs 22 and 23 are fed from the 
outputs 28 and 29 of the circuits 18 and 19 respec- 

so tively. The source 27, and the outputs 28, 29 and 
30 of circuits 18,19 and 20 respectively, are coup- 
led to respective inputs of a voltage summing 
circuit 31 via multipliers 32,33,34 and 35 respec- 
tively. The output of summing circuit 31 constitutes 

55 the output 6 of the generator circuit 7. If the mul- 
tiplication factors of the multipliers 32,33,34 and 35 
are a,0,7, and 5 respectively, and the output volt- 
age of source 27 is V rel , the output voltage of 
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summing circuit 31, i.e. of generator 7, at a time t 
after reset and subsequent release of the integrator 
circuits 18-20 will be V ref (a + 0t + 7 t 2 + 5t 3 ). The first 
two terms correspond to the conventional linear 
voltage ramp starting from a voltage V ref a and 5 
having a slope V ref £. The terms in 7 and 5 are 
quadratic and cubic correction terms respectively; 
quartic etc. correction terms could be added,. or the 
cubic term could be omitted, if desired. 

The multiplication factors a,/9 j7 , and 5 are gen- 10 
erated by digital-to-anaiog converters 36,37,38 and 
39 respectively and fed to inputs 40,41 ,42 and 43 
of the multipliers 32,33,34 and 35 respectively. 
Converters 36 to 39 are fed with the appropriate 
digital data from respective output ports of a digital 75 
data processing arrangement 44 which is consti- 
tuted by a suitably programmed microcomputer. A 
data input port 45 of processor 44 is fed from the 
output 46 of FFT-calculating circuit 15, which out- 
put may also be connected to the data input of a 20 
conventional radar data processing arrangement 
(not shown), as indicated by an arrow 47. Alter- 
natively the conventional radar data processing 
may be performed by the - arrangement 44. The 
reset inputs 24 to 26 of integrator circuits 18 to 20 25 
are fed from an output 48 of arrangement 44. The 
clock inputs 49 and 50 of A/D-converter 11 and 
FFT-caiculator 15 are fed from an output 51 of 
arrangement 44. 

Arrangement 44 is programmed to momentarily 30 
reset integrator circuits 18 to 20 at intervals of 
1mS, except for during the first 10m S of every 
second when it is programmed to maintain these 
circuits continuously in the reset state. Moreover it 
is programmed to apply 2N clock pulses to both 35 
A/D converter 11 and FFT-calculating circuit 15 
during each millisecond, where N is the number of 
points calculated by circuit 15. Thus a 2N-point 
FFT is calculated by circuit. 15 during each mil- 
lisecond and is applied to the input port 45 of aq 
processor arrangement 44. 

Arrangement 44 is programmed to write nomi- 
nal values of a,j3, 7j and 5 to the converters 
36,37,38 and 39 respectively during the initialisa- 
tion phase, where the nominal 7 and 5 may be zero 45 
and the nominal a and 0 are chosen, taking into 
account the output voltage V re i of source 27, the 
time constant of integrator 18, and the frequency 
versus control voltage characteristic of oscillator 1 , 
to give the required starting and finishing frequen- so 
cies respectively of each 1ms frequency sweep 
occurring at the output of oscillator 1 during the 
last 990mS of each second in response to the 
reset signals applied to integrator circuit 18. (It is 
assumed that any value written to any of the con- 55 
verters 36 to 39 is maintained until it is subse- 
quently updated). The starting frequency may lie, 
for example, in the region of 80GHz and the finish- 



ing frequency may be, for example 300MHz higher 
than the starting frequency. 

As mentioned above, the integrator circuits 18 
to 20 are maintained in the reset state during the 
first 10mS of each second, so that the output 
frequency of oscillator 1 is maintained constant at 
the starting frequency during these intervals. Dur- 
ing each of these intervals arrangement 44 is pro- 
grammed to process the ten FFTs then produced 
by circuit 15 to ascertain what the current average 
(noise) signal level is at each of the N outputs of 
circuit 15, and update a record, kept within the 
arrangement 44, of the respective running averages 
of each of these N noise levels in accordance with 
the result. During each 1mS of the remaining 
990mS of each second the arrangement 44 is 
programmed to perform the sequence of oper- 
ations shown in the flow diagram of Figure 2 of the 
drawings, for which purpose it will be assumed, for 
example, that FFT-calculating circuit 15 calculates 
an N = 256-point FFT each time, giving an optimum 
range resolution for reflecting targets of about 0.5 
metre and a maximum range of about 128 metres 
with the sweep rate of 300MHz per mS previously 
quoted for the output signal of oscillator 1 . 

In Figure 2 the various blocks have the follow- 
ing significances 

52 - Start 

53 - Determine that one of the output points 

150 to 256 of FFT-calculating circuit 15 
(corresponding to reflecting targets 
within the range 75 to 128 metres) 
which has the largest amplitude relative 
to the current average noise level for 
this point. 

54 - In the amplitude of this point n at least 

25dB above the current average noise 
level for this point? 

55 - Finish 

56 - Are the amplitudes of points n-3, n-4, n- 

5, n + 3, n + 4 and n+5 relative to the 
respective current average noise levels 
for these points ail at least 3dB less 
than the relative amplitude of point n? 

57 - Determine as a new point n that one of 

the output points 150 to 256 of FFT- 
calculating circuit 15 which has the larg- 
est amplitude relative to the current 
average noise level for this point after 
any point which has already been test- 
ed in a test 54- during the present rou- 
tine, together with its surrounding six- 
teen points, have been excluded from 
consideration. 

58 - Are the amplitudes of points n-6, n-7, n- 

8, n + 6, n+7 and n + 8 all less, than 
10dB above the current noise levels for 
the corresponding points? 
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59- Calculate a 64-point inverse Fast 
Fourier Transform on data the first 
thirty-two points and the final sixteen 
points of which are zero and points 
thirty-three to forty-eight of which are 
the amplitudes of output points n-7 to 
n + 8 respective of the FFT-calculating 
circuit 15 relative to the current average 
noise levels for these points. 

60 - Determine the zero-crossing points of 

the inverse FFT calculated in step 59, 
i.e. positions on the inverse FFT where 
adjacent points have opposite signs, by 
straight-line interpolation between these 
points. 

61 - Determine the maximum absolute value 

of each set of points of the inverse FFT 
calculated in step 59 which lie between 
an adjacent pair of zero-crossing points 
determined in step 60. 

62 - Are the absolute values determined in 

step 61 all within 10% of each other? 

63 - Add a notional cumulative phase 0 ver- 

sus elapsed time t characteristic cor- 
responding to a frequency f 0 to the cu- 
mulative phase 0 versus elapsed time t 
characteristic determined in step 60, 
where f 0 is the frequency corresponding 
to the point n-7 processed in step 59. 

64 - Normalise the characteristic obtained in 

step 63, i.e. divide each of the values of 
cumulative phase 0 by the maximum 
value Omax of cumulative phase, and 
each of the values of elapsed time t by 
the maximum value W of elapsed 
time. 

65 - Fit a cubic polynomial to the character- 

istic obtained in step 64 expressed as 
t/tmax as a function of 0/0 m ax- 

66 - Calculate values of 0,7 and 6 from the^ 

coefficients of the polynomial found in 
step 65, and update respective running 
averages of these values. 

67 - End. 

The processing arrangement 44 is moreover, 
programmed to write the current values of the 
running averages of 0,7 and 5 to the converters 
37,38 and 39 of Figure 1 at intervals of, for exam- 
ple, one minute. 

Thus, each time the routine of Figure 2 is 
performed it is ascertained by means of steps 53 
to 58 whether a radar return having characteristics 
which render it basically suitable for revealing nonr 
linearities in the frequency versus time characteris- 
tic of oscillator 1 is being received by the aerial 3. 
More particularly, steps 53 and 54 select a return 
(if present) which has an appropriately large am- 
plitude, step 56 tests whether the corresponding 



target is no more than three range cells wide at the 
3dB point and step 58 tests whether the cor- 
responding target is isolated in an area at least 
sixteen range cells wide which contains no other 
5 target. If these conditions are satisfied step 59 
effectively constructs the time domain waveform of 
a beat signal between the relevant frequency com- 
ponent of the output signal of mixer 8 and a signal 
having a frequency f 0 corresponding to output point 
70 n-7 of the FFT-calculating circuit 15, the padding of 
the sixteen items of actual data with forty-eight 
zeros having a curve-smoothing function and the 
fact that this "zero-filling" is asymmetrical (the six- 
teen items of actual data being positioned at the 
75 bottom of the upper half of the inversely-trans- 
formed data items) ensuring that the last-men- 
tioned beat signal has a positive frequency. The 
determination of the zero-crossing instants of this 
beat signal (the relevant frequency component of 
20 the output signal of mixer 8 translated to zero 
frequency due to the windowing of sixteen of the 
output points of FFT calculating circuit 11) in step 
60 then gives the times at which successive incre- 
ments of 7T occur in its phase, i.e. gives points on 
25 its cumulative phase 0 against elapsed time t char- 
acteristic. Steps 61 and 62 constitute a check on 
whether the depth of any amplitude modulation of 
the relevant frequency component translated to 
zero frequency is below a predetermined threshold 
30 (because excessive amplitude modulation would 
have produced errors in the positions of the zero 
crossings determined in step 60) and if all is well in 
this respect step 63 translates the cumulative 
phase versus elapsed time characteristic of the 
35 zero-translated frequency component obtained in 
step 60 into a cumulative phase versus elapsed 
time characteristic of the actual relevant frequency 
component of the output signal of mixer 8; a value 
27rf 0 t, where t is the corresponding value of the 
40 elapsed time, is added to each cumulative phase 
value of the first-mentioned characteristic. After 
normalisation of both axes of the resulting char- 
acteristic in step 64 a cubic polynomial 

45 t/W = jS'(0/0max) + YW0max) 2 + *«(0/0max) 3 

is fitted to the resulting set of points (where it will 
be noted that the axes of the characteristic have 
effectively been transposed; the characteristic has 
so effectively been inverted so that it takes the form of 
the corrections required rather than of the errors). 
An appropriate polynomial-fitting algorithm is ob- 
tainable, for example, from the National Algorithm 
Group, (NAG) of Oxford University under reference 
55 E02ACF. This particular algorithm fits the polyno- 
mial in such a way that the worst case deviation of 
any point therefrom is minimised, although of 
course some other criterion for best fit, e.g. 
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minimising the mean square deviation of the var- 
ious points therefrom could alternatively be em- 
ployed, if desired. The result is values for fi\ 7 ' 
and 5'. 

It can be shown that, in order to linearise the 
frequency versus time characteristic of the output 
signal of oscillator 1 in Figure 1 , its normalised 
control voltage (V-V 0 )/(V max -V o ) should vary in the 
same way with normalised elapsed time Vt^ as 
normalised elapsed time t/t max varies with normalis- 
- ed cumulative phase 0/0 max of the selected fre- 
quency component of the output signal of mixer 8. 
(V 0 is the value of the control voltage at the start of 
each frequency sweep, i.e. V ref a). Thus it is re- 
quired in the present case that 

<t/W) 3 
i.e. that 

v-v 0 * a'(v mox -v 0 ) + w'owvoyu, + t 2 7 ^(v mox -v 0 )- 

/(tmax) 2 + t35'{V max -V 0 )/(t max )3 

Thus the coefficients 0, 7 and 5 which it is required 
be applied as multiplying factors to the multipliers 
33,34 and 35 of Rgure 1 are given by 

fi ■ ^'(Vmax-VoVCW) 
7 = 7 , (Vmax-V 0 )/(t max )2 
* = a^max-Vo)/^. 

These coefficients are calculated in step 66 
from the values found for 0\ y and 5' together with 
the values of (V max -V 0 ) and U*. w being pre- 
determined by the duration of each frequency 
sweep of the output signal of oscillator 1 and (V max - 
V 0 ) being predetermined by the required terminat- 
ing frequency of each frequency sweep. 

It was assumed above for the purposes of the 
description with reference to Rgure 2, in particular 
that with reference to steps 64 and 65 thereof, that 
the voltage ramps supplied by the control signal 
generator 7 of Figure 1 to the frequency control 
signal input 5 of radio-frequency source 1 prior to 
each updating of the weighting coefficients applied 
to multipliers 33-35 are each sufficiently linear that 
they can be treated as such for the purpose of the 
calculation of the new values of 0 l7 and 5. (Such an 
assumption enables the amount of data processing 
required to be reduced below that which would 
otherwise be necessary). While such an assump- 
tion is obviously justified initially if, as postulated, 
the values of 7 and 8 are initially set to zero, it may 
not be justified subsequently, especially if the out- 
put frequency versus frequency control voltage 
characteristic of r.f. signal source 1 is itself more 
than slightly non-linear. If this is the case a more 



rigorous calculation of each new set of values for 
0,7 and 5 may be performed in a manner which will 
now be described. 

As mentioned previously, linearization of the 

5 frequency sweeps produced by r source 1 entails in 
general that, for a given frequency sweep from 
source 1, the control voltage V applied to its input 
5 be made to vary in substantially the same way 
with normalized elapsed time t/tmax as it does with 

10 normalized cumulative phase 0/0 max of a beat fre- 
quency signal appearing in operation at the output 
9 of mixer 8. If, for example, correction is required 
up to and including third order as in the example 
described previously, then the above requirements 

is can be expressed mathematically as 

V = «' + ^(t/W) + tWW) 2 + ^(tW (1) 
where a\$\ y % and 8' also satisfy 

20 

V = a'+y S '(0/0 ma x)+ 7 , (0/0mox) 2 +^(0/0 max )3 ( 2 ) 

(If desired Ua* and 0 max may be replaced by t, and 
0 r respectively where 0 r is the value of the cumula- 
25 tive phase at an arbitrary elapsed time tr after the 
beginning of a given control voitage/r.f. output- 
frequency sweep). As mentioned previously, the 
output signal of the control signal generator 7 of 
Rgure 1 is given by 

30 

V=V re( (a + 0t+ 7 t 2 +5t 3 ) (3) 

so in general, for correction up to and including 
third order, it is required that 

35 

V ref a = a' 
Vref£ ~ /37t max 
VrefT = T'/t 2 ™ 

40 and 

VrefS = 5Vt 3 max (4) 

where a\ff, y ' and 5' satisfy expression (2). Thus 
45 the required values of 0, 7 and 5 can be calculated 
by (a) determining the cumulative normalized 
phase 0/0 max versus elapsed time t characteristic of 
a beat signal appearing at the output 9 of mixer 8 
of Rgure 1 during a given frequency sweep of the 
50 output signal of the- r.f. source 1, (b) substituting, 
using equation (3) and the present values of a,0, 7 
and 5, the current corresponding value of V for 
each value of elapsed time t in the characteristic 
determined in (a) to give the corresponding cu- 
ss mulative normalized phase versus control voltage 
characteristic, i.e. pairs of values (V, 0/0 max ), (c) 
fitting a third-order polynomial to these values to 
obtain the coefficients 0V and 8' in equation (2), 
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and (d) using the equations (4) to derive the re- 
quired value of 0,7 and 5 from these coefficients. 
Thus if a more rigorous determination of the re- 
quired values of 0,7 and 5 is required, steps 64 and 
65 of Figure 2 as previously defined may be re- 
placed by 

64 - (a) Normalize the cumulative phases in 

the characteristic obtained in step 63, 
i.e. divide each of the values of cumula- 
tive phase 0 by the maximum value of 
cumulative phase, and (b) substitute 
corresponding values of V for the times 
obtained in step 63. 

65 - Fit a cubic polynomial to the character- 

istic obtained in step 64 expressed as V 
as a function of 0/0 m ax- 
It will be noted that the major data processing 
burden added by the substitution of the above 
operations for the operations 64 and 65 originally 
described with reference to Figure 2 lies in the 
calculation of the corresponding value of V for each 
value of time t in the characteristic obtained in step 
63. If it were assumed that the voltage ramps 
generated by control voltage generator 7 is linear 
prior to the normalization process, then of course 
normalized control voltage (V-Voy^^-Vo) could 
be substituted directly for normalized elapsed time, 
which is the. reason that elapsed time is normalized 
in the step 64 originally described and that this 
step contains no actual substitution of control volt- 
age values for time values. 

As mentioned previously the beat frequency 
which results from reflection by a target at a given 
range will be Doppler-shifted if there is relative 
motion between the target and the radar apparatus. 
There is an implicit assumption in step 63 of Figure 
2 that the frequency f 0 employed therein is not 
significantly. Doppler-shifted. If significant Doppler 
shift is in fact possible its effects can be ascer- 
tained, enabling correction to be made therefor, by 
the known technique of arranging that oscillator 1 
of Figure 1 generates both upwardly-directed and 
downwardly-directed frequency sweeps, these giv- 
ing rise to oppositely directed Doppler shifts in the 
beat frequency due to a given target, and thereby 
enabling the true beat frequency to be ascertained 
by an averaging process. 

It will be evident that many modifications are 
possible to the embodiment described, within the 
scope of the invention as defined by the claims. 
For example, the fact that correction terms up to 
and including third-order are employed is not es- 
sential; even higher-order terms may be employed, 
if desired, or the third order correction may be 
omitted. Obviously the number of integrators and 
associated multipliers included in the control volt- 
age generator 7 of Figure 1 should be chosen in 
accordance with the order of correction provided 



for in each particular case. Moreover, although the 
multipliers 33,34 and 35 are shown being fed from 
respective taps on the cascade arrangement of 
integrators 18,19,20 it will be evident that these 
5 multipliers could alternatively be fed from com- 
pletely separate sets of one, two and three integra- 
tors respectively, which sets could, if desired, be 
fed with input voltages which differ from each other 
and from the voltage fed to multiplier 32. If such 
to different input voltages are employed obviously 
equations (3) and (4) and the calculations employ- 
ing these equations will have to be modified to take 
this into account, i.e. to take into account the 
different values of V re , associated with the respec- 
75 tive weighting coefficients a,fi,y and 5. Furthermore 
the weighting of the various components of the 
control voltage generator output signal may then be 
achieved by adjusting the various V rei if desired, 
rather than by means of the multipliers 33,34 and 
20 35 which may then be omitted. As another exam- 
ple, the screening of the selected frequency com- 
ponent for amplitude modulation (step 62 in Figure 
2) is, although highly desirable, not essential. It 
should be borne in mind though that high degrees 
25 of amplitude modulation (due e.g. to the presence 
of multiple targets at the relevant range) can result 
in large errors when the phase versus frequency 
characteristic is determined in step 60. Indeed, if 
the modulation depth is greater than 100% spur- 
30 ious actual phase reversals will occur. As yet an- 
other example it is not essential to employ asym- 
metrical "zero fill" in step 59, or even "zero fill" at 
all. However, as pointed out previously, "zero fill" 
has a curve-smoothing effect and, if it is made 
35 asymmetrical, it enables positive and negative fre- 
quencies of the zero-translated selected frequency 
component to be differentiated between. An alter- 
native way of achieving such differentiation would 
be, of course, to calculate complex rather than 
40 simple forward and inverse Fourier Transforms. 
However this would necessitate the replacement of 
the single mixer 8 of Figure 1 by a pair of 
quadrature-related mixers and associated output 
channels, so that the asymmetrical zero-fill tech- 
45 nique will often be preferred. (As an alternative to 
positioning the sixteen actual data points at the 
bottom of the upper half of the 64-point transform 
performed in step 59 they could be positioned at 
the top of the bottom half, in which case the 
so frequency f 0 used in step 63 should correspond to 
point n + 8 processed in step 59, rather than point 
n-7). 

The specific figures quoted for various aspects 
of the embodiment described, e.g. number of 
55 points in each FFT calculated by circuit 1 1 , number 
of points in each inverse FFT calculated by ar- 
rangement 44, duration and frequency of each re- 
ceiver noise determining step carried out in re- 
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sponse to the program of the arrangement 44, 
specific criteria employed in selecting a particular 
frequency component of the output signal of mixer 
8, etc. etc., are of course only examples, and may 
be varied as desired in accordance with the pre- 
vailing circumstances. 

Further, in Figure 1 the aerial 3 and the cir- 
culator 4 may be replaced by separate transmit 
and receive aeriais connected respectively to the 
output 2 of the oscillator 1 and to the first input 10 
of the mixer 8. 

From reading the present disclosure, other 
modifications and variations will be apparent to 
persons skilled in the art. Such modifications and 
variations may involve other features which are 
already known in the context of radar apparatuses 
and data processing therein and which may be 
used instead of or in addition to features already 
described herein. Although claims have been for- 
mulated in this application to particular combina- 
tions of features, it should be understood that the 
scope of the disclosure of the present application 
also includes any novel feature or any novel com- 
bination of features disclosed herein either explic- 
itly or implicitly, whether or not it relates to the 
same invention as presently claimed in any claim 
and whether or not it mitigates any or all of the 
same technical problems as does the present in- 
vention. The applicants hereby give notice that new 
claims may be formulated to such features and/or 
combinations of such features during the prosecu- 
tion of the present application or of any further 
application derived therefrom. 

Claims 

1. A method of reducing non-linearity (if any) 
present in the frequency-versus-time charac- 
teristic of the output signal of a swept-fre- 
quency radar transmitter, which output signal 
is derived from a radio-frequency source the 
output frequency of which is dependent on the 
value of a control signal applied in operation to 
a control signal input of said source, the output 
of a control signal generator being coupled to 
said control signal input for applying said con- 
trol signal thereto to repeatedly sweep the 
frequency of said output signal, in which meth- 
od said output signal is mixed with a delayed 
version thereof to yield a beat signal and said 
characteristic is adjusted in accordance with 
any variation of. a frequency of this beat signal 
occurring, during individual frequency sweeps 
of said output signal, characterized in that said 
control signal is formed by combining a time- 
independent component and n individually 
weighted time-dependent components, where 
n is greater than one, where the first time- 



dependent component is derived from the out- 
put of a single resettable analogue integrator 
circuit which is fed from the output of a. d.c. 
reference source and where the mth time-de- 

5 pendent component (m = 2 lv ..n) is~derived from 

the output of a. cascade combination of m 
resettable analogue integrator circuits which Is 
fed from the output of a d.c. reference source, 
in that the delayed version of the transmitter 

w output signal is constituted by a version of this 

output signal which has been transmitted and 
reflected back by a target or targets, in that the 
cumulative-and-normalized-phase versus time 
characteristic of a specific spectral component 

75 of said beat signal is determined for an individ- 

ual frequency sweep of said transmitter output 

signal, in that the coefficients of the first nth 

order terms of an nth order polynomial which 
fits the dependence of said control signal on 

20 said cumulative and normalized phase are de- 

termined from said cumulative-and-normalized- 
phase versus time characteristic and the con- 
comitant variation of said control signal with 
time, and in that the coefficients so determined 

25 are used to update the weighting coefficients 

of the first....nth said time-dependent compo- 
nents respectively to reduce any said non- 
linearity. . 

30 2. A method as claimed in Claim 1, wherein said 
beat signal is analysed into its spectral compo- 
nents by calculating a Fourier Transform of 
said beat signal, and data defining said spe- 
cific spectral component is derived from said 

35 Fourier Transform by calculating an inverse 

Fourier Transform of a windowed portion there- 
of. 

3. A method as claimed in Claim 2, wherein the 
40 calculated inverse Fourier Transform is of said 

windowed portion supplemented with zero 
data. 

4. A method as claimed in Claim 3, wherein the 
45 supplementation with zero data is asymmet- 
rical. 

5. A method as claimed in Claim 2, wherein said 
windowed portion is selected on the basis that 

so a single isolated frequency component is in- 
dicated therein. 

6. A method as claimed in Claim 3, wherein said 
windowed portion is selected on the basis that 

55 a. single isolated frequency component is in- 

dicated therein. 

7. A method as claimed in Claim 4, wherein said 
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windowed portion is selected on the basis that 
a single isolated frequency component is in- 
dicated therein. 

Radar apparatus comprising a radio-frequency 
source having an output coupled to signal op- 
erating means for radiating a transitter signal 
and for receiving at least a return signal, a 
control signal generator having an output coup- 
led to a frequency control signal input of said 
radio-frequency source for applying a frequen- 
cy control signal to said input to repeatedly 
sweep the frequency of the output signal of 
said source, and a mixer having a first input 
coupled to the signal operating means for re- 
ceiving a version of the output signal of said 
source which has been transmitted and re- 
flected back by a target or targets, a second 
input coupled to the output of said radio-fre- 
quency source for receiving a sample of the 
output signal of said source as a reference 
signal, and an output for a beat signal between 
said version of the output signal of said source 
and said sample, said control signal generator 
comprising component combining means for 
combining a time-independent component and 
n individually weighted time-dependent com- 
ponents to form said frequency control signal, 
where n is greater than one, a resettable ana- 
logue integrator circuit arrangement for gen- 
erating said time-dependent components, and 
means for individually weighting each said 
time-dependent component, said component 
combining means having a respective input 
corresponding to each said time-dependent 
component and said resettable analogue in- 
tegrator circuit arrangement comprising a sin- 
gle resettable analogue integrator circuit in- 
cluded between the output of a d.c. reference 
source and the combining means input for the 
first time-dependent component and a cascade 
combination of m resettable analogue integra- 
tor circuits included between the output of a 
d.c. reference source and the combining 
means input for the mth time-dependent com- 
ponent (m = 2 n), Ihe apparatus further in- 
cluding means for determining the cumuiative- 
and-normalized-phase versus time characteris- 
tic of a specific spectral component of the 
mixer output signal for an individual frequency 
sweep of the radio-frequency source, means 
for determining the coefficients of the 
first nth order terms of an nth order polyno- 
mial which fits the dependence of said fre- 
quency control signal on said cumulative-and- 
normalized-phase from said cumulative-and- 
normalized-phaseversus time characteristic 
and the concomitant variation of said control 



signal with time, and means for updating the 

weighting coefficients of the first ,nth said 

time-dependent components in accordance 

with the coefficients of the first nth order 

5 terms respectively so determined to reduce 

non-linearity (if any) present in the frequency 
sweeps of the output signal of said radio- 
frequency source. 

10 9. An apparatus as claimed in Claim 8, further 
comprising analysing means for analysing said 
beat signal into its spectral components by 
calculating a Fourier Transform of said beat 
signal and means for calculating an inverse 

75 Fourier Transform of a windowed portion of the 

Fourier Transform to derive data defining said 
specific spectral component. 
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